Ferromagnetic materials in the zinc-blende structure 
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New materials are currently sought for use in spintronics applications. Ferromagnetic materials 
with half metallic properties are valuable in this respect. Here we present the electronic structure 
and magnetic properties of binary compounds consisting of 3d transition metals and group V ele- 
ments viz. P, Sb and As in the zinc-blende structure. We demonstrate that compounds of V, Cr and 
Mn show half metallic behavior for appropriate lattice constants. By comparing the total energies 
in the ferromagnetic and antiferromagnetic structures, we have ascertained that the ferromagnetic 
phase is stable over the antiferromagnetic one. Of the different compounds studied, the Cr based 
systems exhibit the strongest interatomic exchange interactions, and are hence predicted to have 
the highest critical temperatures. Also, we predict that VAs under certain growth conditions should 
be a semiconducting ferromagnet. Moreover, critical temperatures of selected half metallic com- 
pounds have been estimated from mean field theory and Monte Carlo simulations using parameters 
obtained from a ab-initio non-coUinear, tight binding linearized mufBn-tin orbital method. ^From a 
simple model, we calculate the reflectance from an ideal MnAs/InAs interface considering the band 
structures of MnAs and InAs. Finally we present results on the relative stabilities of MnAs and 
CrSb compounds in the NiAs and zinc-blende structures, and suggest a parameter space in substrate 
lattice spacings for when the zinc-blende structure is expected to be stable. 

PACS numbers: 71.20.Be,75.50.Pp,71.70.Gm 



I. INTRODUCTION 



For magnetoelectronic devices, highly spin polarized 
ferromagnets have a large potential in applications 0. 
Half metallic ferromagnets (HMF) are the most desirable 
materials in this respect. These materials have one spin 
channel that is metallic, while the other spin channel has 
a gap at the Fermi level. Hence, 100 % spin polarization 
of the conducting electrons is expected from these ma- 
terials. Based on band-structure calculations, de Groot 
et al. first predicted half metallic behavior of Clb 
type Heusler alloys of NiMnSb and PtMnSb. Since then, 
there have been predictions of many half metallic mate- 
rials. A practical aspect of these materials is the recent 
breakthrough utilizing the spin of electrons in existing 
semiconductor electronic devices. To manipulate spins 
of electrons inside the semiconductors, one first needs to 
inject spins inside the semiconductors. Spin injection can 
be done either from ferromagnetic metals or from diluted 
magnetic semiconductors (DMS) doped with magnetic 
ions. Ohno et al. 3j successfully injected spin polarized 
holes from a p-type HI-V DMS to a semiconductor. From 
theoretical calculations of magnetotransport, Egues Q 
proposed a spin filter for spin-dependent optoelectron- 
ics applications. Later on, Fiederling et al. Q success- 
fully injected spin polarized electrons from an n-type II- 
VI semiconductor. However, according to Schmidt et al. 
0, the spin injection from a metal to a semiconductor is 
practically impossible due to the huge conductivity mis- 
match between a metal and a semiconductor. The use of 
a DMS is instead rather advantageous. Mn doped GaAs 
is a prototype of this class of systems but the highest 
available Curie temperature is only 110 K Q due to the 
limited solubility of Mn in GaAs. The use of other half 



metallic systems e.g. Cr02 and Lao.rSro.sMnOs is also 
restricted due to lower Curie temperature. Very recently, 
there have been reports on room temperature ferromag- 
netism in Mn doped GaP compounds §\. Hence, the 
quest for new half metallic materials with room temper- 
ature ferromagnetism continues. 

Recently, there have been efforts to grow binary com- 
pounds in the zinc-blende structure using molecular 
beam epitaxy. One of the advantages of these materials is 
that the zinc-blende structure is maintained when grown 
on a zinc-blende semiconductor substrate. Akinaga et al. 
grew CrAs on GaAs substrates and found CrAs to be a 
ferromagnet having a Curie temperature higher than 400 
K. The experimental results were complemented by ah- 
initio electronic structure calculations showing that the 
magnetic moment was 3 ur. Also CrSb was grown in 
the zinc-blende structure [iQ- These binary compounds 
do not exist naturally in the zinc-blende structure but 
the non equilibrium growth by MBE has facilitated the 
creation of these compounds with novel magnetic prop- 
erties. 

Sanvito and Hill [Til l examined the ground state prop- 
erties of MnAs in a NiAs and a zinc-blende structure s by 
first principles pseudopotential calculations. Shirai |l2| 
performed electronic structure calculations for the zinc- 
blende compounds of transition metal (TM) and As. The 
early TM-As compounds were shown to be ferromagnets 
and the later compounds were antiferromagnets. Also, 
electronic structure calculations have been performed for 
different crystal structures of MnP, MnAs and MnSb by 
Continenza et al. Galanakis |0| examined the half 
metallicity of the bulk and surface of zinc-blende CrAs in 
various lattice constants by electronic structure calcula- 
tions. A full potential screened Korringa-Kohn-Rostoker 
method was used. Galanakis concluded that a Cr ter- 
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minated CrAs (001) surface remains half metallic both 
for the experimental lattice constants of GaAs and InAs. 
On the contrary, an As terminated surface loses its half 
metallicity due to the presence of the surface dangling 
bond states. 

The motivation of this work is to design new materials 
in a crystal structure that is suitable for spin electron- 
ics applications and provide a guideline for experimental 
work that attempt to grow these materials. The paper 
is organized as follows : In section 'Computational de- 
tails', we describe in detail the computational part. In 
the results section, a subsection 'Electronic structure and 
magnetism' describes the electronic structure and mag- 
netism of the binary compounds considered followed by 
the discussions about the interatomic exchange interac- 
tions and critical temperatures in some of the systems 
studied. Finally, structural stabilities of MnAs and CrSb 
compounds in various structural phases are presented. 



II. COMPUTATIONAL DETAILS 

A zinc-blende structure has two fee sublattices dis- 
placed from each other along the body diagonal by a vec- 
tor (1:11:). For calculations in the ferromagnetic phase, 
the unit cell used in the calculations is the same as in the 
fee lattice and has 1 TM at the (000) and the group V 
atom at the (^i^:) position. Calculations have been per- 
formed by an ab-initio plane wave pseudopotential code 
(VASP) [13 . Pseudopotentials generated using the Pro- 
jector augmented wave (PAW) 16] method were used 
throughout the calculations. These pseudopotentials are 
suitable for transition metals (TM) and the calculations 
yield accurate results |3 . A kinetic energy cut-off of 600 
eV was used for the plane waves included in the basis 
set. Perdew-Wang GGA (generalized gradient approxi- 
mation) exchange -correlation was used. According 
to the calculations of Continenza et at T^, GGA is 
essential for obtaining accurate equilibrium volume and 
magnetic moments. A 20x20x20 k-points grid was used 
in the Monkhorst Pack scheme fig^ which yielded 770 k 
points in the irreducible Brillouin zone and subsequently 
used for the Brillouin zone integration within the tetra- 
hedron method. Local properties such as local density 
of states and local magnetic moments were calculated by 
projecting the wave functions into spherical harmonics as 
described in ref. |2Q|| . The values of the experimental lat- 
tice constants used in these calculations were taken from 
Ashcroft and Mermin . Calculations were done for the 
TM-V compounds for the lattice constants of the III-V 
semiconductors having the same Group V element. In 
summary, we considered lattice constants of GaAs, GaP, 
GaSb, InAs, InP and InSb compounds which covers a 
wide range (5.45 A for GaP to 6.48 A for InSb). 

The pseudopotential calculations have been comple- 
mented [231 with calculations based on non-coUinear 
tight-binding linear mufhn-tin orbital method (TB- 
LMTO- AS A) j2^ including combined correction terms in 



the one-electron Hamiltonian'24'|. The generalized gradi- 
ent approximation (GGA) were used for the exchange- 
correlation potential 25] . Equal Wigner-Seitz sphere 
radii were used for all atoms, as well as for the empty 
spheres representing the two types of interstitial sites 
in the zinc-blende structure. A dense 27x27x27 k-point 
grid was used in the Brillouin zone integration and to 
speed up convergence each energy eigenvalue was asso- 
ciated with a Gaussian of width 85 meV. Spin wave 
excitations have been studied using the frozen-magnon 
approachlH IH Then the total energy £'(q,6') is 
calculated for spiral magnetic structures where each mag- 
netic moment is defined by the Euler angles 



const; (f)i = q ■ Tii 
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where 9i and (f>i are the polar and azimuthal angles, 
the position of ion i and q is the spiral propagation 
vector. S pin spirals possess a generalized translational 
symmetry |29|, which allow for calculations in the chemi- 
cal unit cell without need for a supercell. It can be shown, 
that the true spin wave excitation energies are related to 
the spiral energies through the relation 



4 AE{q,t 
M sin? 9 



(2) 



where A£;(q, 9) = £'(q, 9) - E{0, 9). The last equation is 
valid for one atom per cell. Thus the formula is correct 
if one consider the magnetic response from the transi- 
tion metal only and neglect the small induced magnetic 
moment on the group V element. In our calculations, 
we use an angle 9 — 20° and the spiral energies were 
obtained by employing the magnetic force theorem. We 
tested other angles as well as self-consistent calculations 
and found that the differences in total energy is within a 
few percent. 

In order to obtain thermodynamic properties, the spin 
wave excitation energies were mapped to an effective 
Heisenberg Hamiltonian of classical spins 



(3) 



where are exchange interactions between the tran- 
sition metal atoms. The critical temperatures Tc were 
estimated both from mean field theory^ (MFA) and 
Monte Carlo simulations (MC). For the MC simulations, 
we used the standard single flip Metropolis algorithm and 
the critical temperature were obtained from the 'cumu- 
lant crossing method' [s^]. In this method, the reduced 
fourth order cumulant of the order parameter (magneti- 
zation) U L , defined as 



(4) 
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is calculated for different sizes of the simulation box. The 
curves oiU^ have a common intersection point at a fixed 
point [/*, which corresponds to the critical point. Hence, 
a value of Tc is obtained from the intersection point of 
Ul for different lattice sizes. 



III. RESULTS AND DISCUSSIONS 
A. Electronic structure and magnetism 

The zinc-blende lattice introduces a tetrahedral crys- 
tal field under which the transition metal atomic d-states 
are split into triply degenerate t2 and doubly degener- 
ate e states, e states lie lower in energy than t2 states 
in a tetrahedral crystal field in contrast to the octahe- 
dral crystal field. The t2 states hybridize with the anion 
p-states. As a result, bonding and antibonding states 
appear. Also, depending on the exchange splitting, the 
bonding and antibonding states may lie in different po- 
sitions for different spin channels, e states are of non- 
bonding character and do not take part in the hybridiza- 
tion. 

In Fig. l(a-f), we present the spin resolved DOSs of 
V, Cr and Mn compounds having the lattice constants of 
GaAs, InAs, GaSb, InSb, GaP and InP. In figs, (a), (b), 
(c), (d) and (f), all the compounds show half metallic 
behavior since there is no state at the Fermi level for the 
spin down DOS. The only exception is the case depicted 
in fig. 1(e) which doesn't exhibit half metallic ferromag- 
nctism as the Fermi level cuts both the spin bands. For 
the lower lattice constant of GaP, the p — d overlap and 
hybridization of the spin down states broadens the band 
to cut the Fermi level. Also, as the width of the gap in 
the spin down DOS is quite large, the Fermi level stays 
within the gap with the variation in the lattice constant. 
It is worthy to be mentioned here that these compounds 
can be regarded to be the derivatives of diluted magnetic 
semiconductors, e.g. Mna;Gai_2;As for a Mn doped GaAs 
system having a;=l. Several first principles calculations 
[31I Is^ Is^ 1^ confirmed that the total spin magnetic 
moment of Mn doped GaAs for various concentrations of 
Mn is 4 hb I cell with a half metallic nature. It was also 
shown [13 that Cr doped GaAs is a half metallic ferro- 
magnet with a total spin magnetic moment of 3 /is /cell. 
We find here a similar trend for the extreme limit at a 
high TM concentration. A systematic variation in the 
DOSs of V, Cr and Mn is seen. As the number of valence 
electrons increases by going from V to Mn, the main fea- 
tures of the DOSs shift to lower energies relative to the 
Fermi level. This can be seen from figs. 1 (b) and 1 (f). 
The exchange splitting is maximum for Mn and here one 
finds a maximum separation between the spin up and 
spin down bands. 

A close inspection of Fig. 1(a) reveals that VAs be- 
haves like a magnetic semiconductor for the lattice con- 
stant of GaAs, i.e. the Fermi level passes through a gap 
for both spin up and spin down channels. We varied 



_ (a) M 


VAs 

CrAs 


""Ml 


VAs 

CrAs 

— MnAs 


- (GaAs) 

1 , 1 




i 

1 




(InAs) 


' 1 ' 

(c) 


\ 

it 


1 \ , \ , 

VSb 

l^^^^^^^^-^^^CrSb 


1 ' Si 
(d) _ 11 1 


— 1 1 1 1 — 

VSb 

CrSb 


- (GaSb) 1 




■ ! 


(InSb) 


— 1 — 1 — 1 — 1 — 


— 1 — 1 — 1 — 1 — 

VP 

CrP 


— , 1 ! — i — 
_(f) \ '\ ill 


— 1 1 1 1 — 

VP 

CrP 
MnP 


(GaP) \ 










(InP) 



-5 -3-113-5-3-1135 
E - (eV) 

FIG. 1: Spin resolved total (sum over atoms and orbitals) 
density of states (DOS) calculated for the unit cell. In (a), 
we show the DOSs for VAs and CrAs in the GaAs lattice 
constant, in (b), VAs, CrAs and MnAs in the InAs lattice 
constant, in (c), VSb and CrSb in the GaSb lattice constant, 
in (d), VSb and CrSb in the InSb lattice constant, in (e), VP 
and CrP in the GaP lattice constant and in (f), VP, CrP and 
MnP in the InP lattice constant. The energies are shifted 
with respect to the corresponding Fermi levels and energies 
are in electron volts. 

the lattice constant of VAs around the GaAs lattice con- 
stant and found that the small gap is sensitive to the 
lattice spacing. Only in a very small region around the 
GaAs lattice constant, do we find that the gap persists. 
This is an interesting result since it is one of quite few 
ferromagnetic semiconductors found in a stoichiometric 
composition. In addition, it is to be noted that unlike 
CrAs, MnAs doesn't show half metallic behavior for the 
GaAs lattice constant. It is instead half metallic for the 
InAs lattice constant (see inset of Fig. 6(b)). 

In Table I, we present the total magnetic moment/cell 
for all the compounds considered. The half metallic fer- 
romagnetic moments are marked in bold. Half metallic- 
ity may be monitored from the integer moments of the 
unit cell. It is seen that from all the combinations, VAs 
and CrAs are HMFs for both GaAs and InAs lattice con- 
stants. VSb and CrSb are also half metallic for GaSb and 
InSb lattice constants. VP, CrP and MnP are only HMFs 
at the InP lattice constant. This shows that the lattice 
constants and the resulting electronic structures play an 
important role in determining half metallicity, a finding 
that is also evident from Fig. 1. The moments listed 
in Table I are consistent with the calculated moments of 
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Cr and Mn impurities in a GaAs host, if one neglects 
the effect of other defects e.g. As antisites [s^]. To be 
more specific we note that Mn as an impurity in GaAs 
carries a magnetic moment of 4 //s/Mn atom in all theo- 
retical first principles calculations based on local density 
approximation (LDA) that neglect other lattice defects 
such as As antisites. If one considers As antisite de- 
fects, or Mn interstitial atoms that couple antiferromag- 
netically to the Ga substituted Mn atoms, the effective 
moment becomes ^2.4 /i^/Mn atom in agreement with 
experimental data 34] . For the currently investigated 
materials only a limited set of experimental data have 
been published. However, CrAs has been synthesized in 
the zinc-blende structure and the magnetic moment was 
measured to be 3 /is /Cr atom 9] in agreement with the 
present theory and previous calculations 9] . Since Mn 
has one extra electron in the spin up band, it is natural 
that the moment of MnAs is 4/is/Mn atom, which is in 
agreement with the moment of dilute Mn in GaAs (if one 
neglects defects). 

The explanation of the half metallicity and the integer 
moments is similar to that of a III-V diluted magnetic 
semiconductor. In a III-V DMS, a cation is removed 
by creating a vacancy containing 3 holes in the valence 
band. When a transition metal occupies this cation site, 
it donates 3 electrons to the anion dangling bonds to sat- 
urate the bonding. The remaining unpaired d electrons 
of the TMs are responsible for the magnetic moments. 
Due to the large exchange splitting according to Hund's 
rule, the spin up and spin down d bands are separated 
from each other having different degree of hybridization 
with the valence bands. Hence, V, Cr and Mn have mag- 
netic moments of 2, 3 and 4 /i^, respectively from the 
unpaired spin up d electrons. However, the widths of the 
d-bands are fairly large indicating a substantial amount 
of hybridization with the anion-p orbitals. 

TABLE I. Total magnetic moment/cell (in /is). Each 
entry corresponds to the binary alloy consisting of a tran- 
sition metal from the leftmost column and the group V 
element in the subsequent columns. The bulk semicon- 
ductors whose lattice constants are considered for each 
calculation are shown in suffices. 
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We also investigated compounds of Ti, Fe, Co and Ni, 
and found that they do not show half metallic ferromag- 
netism (see Table I). Fe compounds have fairly big mo- 
ments in the ferromagnetic phase except for FeP in the 
GaP lattice constant. CoSb in the InSb lattice constant 
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FIG. 2: Exchange interactions for V, Cr and Mn compounds. 
Note tliat the data are given for lattice constants correspond- 
ing to GaAs, InAs, GaSb, InSb and InP substrates. 



has the highest moment among the Co based compounds. 
Ti and Ni based compounds do not have prominent fcr- 
romagnetism. As only V, Cr and Mn based compounds 
show interesting and useful magnetic behavior, we will 
now concentrate on those compounds only. 

TABLE II. Total energy differences in eV/Mn atom 
{AE = Eafm-Efm) where FM and AFM are the fer- 
romagnetic and antiferromagnetic arrangements between 
Mn atoms respectively. The symbol '-' is for the compo- 
sitions where calculations were not performed. 
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B. Interatomic exchange interactions 

In this section, we compare the total energies of the 
HMFs found from our calculations, in the ferromagnetic 
(FM) and antiferromagnetic (AFM) phases. For brief- 
ness, calculations were only made for selected V, Cr and 
Mn compounds (see Table II). In our calculations, we 
considered an 8 atom (TM)4X4 cell in a layer geome- 
try in the (001) direction with 2 TM (TM/) atoms in 
the plane at z=0.0 and 2 TM atoms at z=a/2 (TM//), 
where a is the lattice parameter. Here TM and X in- 
dicate transition metal (V, Cr and Mn) and the group 
V element (As, P and Sb) respectively. TM/ and TM// 
atoms can couple ferromagnetically or antiferromagnet- 
ically. The energy differences are given in Table II and 
are also plotted in Fig. 2. The energy differences, which 
refiects the interatomic exchange coupling strengths are 
fairly large and are found to always stabilize ferromag- 
netism. From Fig. 2 and Table II, it is clear that the 
strongest exchange interaction is found for the Cr based 
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FIG. 3: Spin wave spectra along high-symmetry lines in the 
Brillouin zone of a) VAs b) CrAs c) MnAs. Circles are cal- 
culated values for the lattice constant of GaAs and squares 
are for the lattice constant of InAs. Solid lines represent the 
fitted data. 

compound, where the energy difference is quite apprecia- 
ble, indicating a high Curie temperature. An inspection 
of the DOSs for the FM and the AFM structures reveal 
that for the AFM structure, the DOS at the Fermi en- 
ergy is non zero for both spin up and spin down bands 
whereas in the FM alignment, only spin up states exist 
at the Fermi level. This stabilizes the FM alignment, 
since the eigenvalue- sum in general lowers the total en- 
ergy when electron states are removed from the Fermi 
energy. 

In this section, we also present results of the calcu- 
lations of critical temperatures using spin wave spectra 
and statistical mechanics. We have calculated Tc for VAs, 
CrAs and MnAs for lattice constants of GaAs and InAs 
by extracting exchange interactions from calculated spin 
wave spectra, which subsequently were used in a classi- 
cal Heisenberg Hamiltonian. In Fig. 3, the calculated 
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magnetic and is therefore not shown in the figure. The 
other systems were found to be ferromagnetic and half 
metallic with total magnetic moments of 2.0, 3.0 and 4.0 
/.ts, respectively for VAs, CrAs and MnAs. In Table. Ill, 
calculated values of Tc are shown. In the case of CrAs, a 
very high Tc around 1000 K is estimated which could be 
very promising from an application point of view. T*^*-^ is 
smaller than T^ '^^ due to the fact that fluctuations are 
absent in the latter. Normally, MC simulation gives Tc 
values closer to experimental values than those obtained 
from mean field theoryjU. 

TABLE III. Calculated critical temperature of VAs, 
CrAs and MnAs from the MFA and MC simulations. 
Values at the left correspond to the lattice constant of 
GaAs and the right correspond to InAs. 
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As is clear from Table III, the critical temperatures of 
VAs and CrAs are in general larger for the substrate with 
the smallest lattice constant (GaAs). Hence it is clear 
that the critical temperatures of VAs and CrAs follow the 
same trend as that of rare-earth systems, where in gen- 
eral a reduced lattice constant results in enhanced criti- 
cal temperatures (as for instance measured by an applied 
external pressure 1^3 )• In rare-earths this phenomenon 
is easily understood from the RKKY interaction, which 
is known to be responsible for the magnetic ordering in 
these elements. By analogy we can say that the presently 
studied materials (VAs and CrAs) have rare-earth like 
characteristics, in the sense that the exchange interaction 
between localized moments is mediated by more diffuse 
valence states. 

In fig. 4, we show the exchange interactions (Jy) in 
CrAs as a function of Cr-Cr distance. It is clear that the 
first nearest neighbor exchange interaction is the domi- 
nant one and is ferromagnetic in nature. Jij decreases 
rapidly with increasing distance between Cr atoms. Fur- 
ther neighbor exchange interactions are antiferromag- 
netic in nature and are very weak. Similar kind of ob- 
servation was reported by Sandratskii and Bruno [36| in 
the case of Mn doped GaAs. 



IV. BAND STRUCTURE AND CALCULATION 
OF REFLECTANCE 



FIG. 4: Exchange interactions between the Cr atoms in CrAs 
as a function of Cr-Cr distance (in units of GaAs lattice con- 
stant). 

spin wave spectra for the systems are shown. MnAs in 
the lattice constant of GaAs was found to be antiferro- 



In Fig. 5(a) and (d), we show the band structures 
of MnAs in the InAs lattice constant for two symmetry 
directions viz. F-L and F-X directions respectively. In 
(b) and (c), we show the same for bulk InAs calculated 
for the experimental lattice parameter. Our calculation 
does not produce a gap at the F point, so the top of 
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the valence band and the bottom of the conduction band 
merge with each other. This is a well known limitation 
of using LDA or GGA within density functional theory. 
For MnAs, the spin up bands cross the Fermi level for 
the two symmetry directions shown the figure and also 
for other directions (omitted in the figure). There is gap 
in the spin down channel. As described above, the e 
bands are dispersionless as they do not mix with others. 
The antibonding t2 bands show large dispersion as they 
are hybridized strongly with the As p bands. 

MnAs InAs InAs MnAs 
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FIG. 5: Band structure of MnAs and InAs in two symmetry 
directions. For MnAs, the filled up triangles and empty down 
triangles indicate spin up and spin down bands respectively. 
The horizontal lines indicate the Fermi level of MnAs and the 
valence band maximum of InAs. 

We now give a rough estimate of the re- 
flectance and transmittance from a ferromagnetic 
metal/semiconductor interface from the band matching 
of the individual species. We use a simple free electron 
model, that was used also by Kilian and Victora [s^ . 
For an electron traveling through a one-dimensional 
potential step, the reflectance is given by 

^ ^ (fc/m - ksc)'^ 

where kfm and ksc are the wave vectors in the ferromag- 
netic metal and semiconductor respectively. As the spin 
up bands are the only states available at the Fermi en- 
ergy, we consider the following for the spin up electrons 
only. For InAs, as the conduction band minimum is at 
the r point, the lowest energy levels available to accept 
spins from the metal are at the F point. We calculate 
R[100], R[UO] and i?[Ill] for [100] (F-A), [110](F-X) 
and [1I1](F-L) directions, respectively. It is known from 
experiments that for e.g. an interface of Au/n-InAs, 
the Fermi level is pinned 0.1 eV above the conduction 
band minimum. We use this information in our model 



assuming that the properties of the metal do not alter 
strongly the band alignments. Also, the spins within 
{ksT 0.025ey) of the MnAs Fermi level are the only 
ones available for transport properties. From this infor- 
mation we can calculate ksc for the three different direc- 
tions, kfm foi' the three different directions are available 
from the band structures of MnAs directly. Then we 
calculated the reflectance from equation 1. The calcu- 
lated values are i?[100]=0.73, i?[110] = 0.89 and i?[lll] 
= 0.71. This means that the transmissions in these di- 
rections will be rather small. However, one should bear 
in mind that in reality the interface is not ideal. There 
arc interface states and different scattering mechanisms 
which can further decrease or increase the transmission 
through the interface. A full transport calculation taking 
into account a proper model of the interface could give a 
more realistic estimate, and our model should be viewed 
as an estimate only. 



V. STRUCTURAL STABILITY 

In this section, we present the relative structural sta- 
bilities of these compounds. For briefness, we choose to 
illustrate our conclusion only for 2 out of the 9 HMFs 
shown in Table I, i.e. CrSb and MnAs. The conclusions 
reached for these two compounds should be valid for all 
the HMFs listed in Table I. We compare the total en- 
ergies of two crystal structures, hexagonal NiAs (that is 
the low temperature structure) and the zinc-blende (ZB) 
structure. 
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j I , 1 , 1 , 1 , 1 -14.25 I ' ' ' ' ' 1 

30 40 50 60 70 25 35 45 55 

Volume/f.u. (Ang.^) Volume/f.u. (Ang.^) 



FIG. 6: Total energy vs. volume for (a) CrSb in the zinc- 
blende and the NiAs structure (b) MnAs in the zinc-blende 
and the NiAs structure with (inset) variation of total magnetic 
moment/cell with lattice constant for zinc-blende MnAs. 

From Fig. 6, it is evident that for both MnAs and 
CrSb, the NiAs structure is more stable than the ZB 
structure. Similar results were obtained for MnAs by 
Sanvito and Hill (TJ- This is in accordance with nat- 
urally occurring compounds of CrSb and MnAs in the 
NiAs structure. However the total energy curves for 
the NiAs structure crosses the total energy curves of ZB 
structures for certain volumes. This occurs for a volume 
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of 46.3 A 3/f u in MnAs and 59.6 A ^/i.u. in CrSb. In- 
cidentally, these volumes are quite close to that of GaAs 
in case of ZB-MnAs and that of GaSb in case of ZB-CrSb 
compounds. For larger volumes, the ZB structure is more 
stable than the NiAs structure. Also, these compounds 
become half metallic ferromagnets according to the calcu- 
lations. Therefore, ZB heterostructures of suitable com- 
binations should be possible to use to obtain half metal- 
lic materials with a large critical temperature. Hopefully 
this fact can be a useful guideline for experimental work 
in trying to grow heterostructures for spin-electronics ap- 
plications. Finally, the variation of magnetic moment of 
MnAs with the variation in lattice constant is shown in 
the inset of Fig. 6(b) where it is clear that it becomes 
half metallic only near the vicinity of a lattice constant 
of 6 A . 

The results presented above are in excellent agreement 
with the existing ones, wherever available. However, 
there are certain limitations of our approach. One of 
them is the underestimation of band gaps in semicon- 
ductors within the use of LDA or GGA. A more appro- 
priate way is to use GW approximation [s^ for these 
materials which corrects the band gaps extremely well. 
To our knowledge, no theoretical calculation exists for 
diluted magnetic semiconductors using this more com- 
plicated theory. The other issue is the effects of native 
defects which are always present in the growth processes. 
One should also take into account these effects to model 
a more realistic situation. But we believe that the re- 
sults presented in this paper are quite significant as far 
as trends are concerned and also a guideline for future 
experiments. 



VI. CONCLUSION 

In conclusion, we have identified a number of possible 
half metallic ferromagnets in the zinc-blende structure. 
Ferromagnetism has been shown to be more stable than 
an antiferromagnetic structure and exchange couplings 
are calculated to be sufficiently strong to produce high 
Curie temperatures, especially for Cr based compounds. 
We also give a parameter range in substrate lattice con- 
stants where MnAs and CrSb are expected to grow and 
to have half metallic properties. Finally, we predict that 
VAs grown at suitable conditions should be a ferromag- 
netic semiconductor with a large critical temperature. 
We hope to motivate experimental work in trying to grow 
such materials for future applications. From a theoretical 
point of view, the future scope of this work is to investi- 
gate the interface properties of these binary compounds 
and the zinc-blende semiconductors e.g. CrAs and GaAs. 
A detailed investigation including structural relaxations 
of these interfaces is in progress. 
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